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Abstract 

We studied the low speed fracture regime (10~ 4 m.s" 1 - 1(T 9 m.s 1 ) in different 
glassy materials (soda-lime glass, glass-ceramics) with variable but controlled length 
scale of heterogeneity. The chosen mechanical system enabled us to work in pure 
mode I (tensile) and at a fixed load on DCDC (double cleavage drilled compression) 
specimen. The internal residual stresses of studied samples were carefully relaxed 
by appropriate thermal treatment. By means of optical and atomic force (AFM) 
microscopy techniques fracture surfaces have been examined. We have shown for the 
first time that the crack front line underwent an out-of-plane oscillating behavior 
as a result of a reproducible sequence of instabilities. The wavelength of such a 
phenomenon is in the micrometer range and its amplitude in the nanometer range. 
These features were observed for different glassy materials providing that a typical 
length scale characterizing internal heterogeneities was lower than a threshold limit 
estimated to few nanometers. This effect is the first clear experimental evidence of 
crack path instabilities in the low speed regime in a uniaxial loading experiment. 
This phenomenon has been interpreted by referring to the stability criterion for a 
straight crack propagation as presented by Adda-Bedia et al. [1]. 
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1 Introduction 



Quasistatic brittle fracture is a subject of importance for both practical and 
fundamental reasons. Two classes of experimental set-up have been primarily 
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developed in order to propagate cracks in a stable and controllable manner. 
The first one is based on a thermally induced stress field [2-4]. An other classi- 
cal way to propagate slow cracks is to use double cleavage drilled compression 
(DCDC) specimen. This method initially developed by Janssen [5] has numer- 
ous advantages: compression loading, mid plane crack stability, auto precrack- 
ing and pure mode I (tensile) crack tip opening. However numerous examples 
of unstable straight crack propagation were evidenced in thermal experiments 
even in the low speed regime. According to values of thermal gradient and 
width of the glass plate, three types of crack propagation regime can develop 
from the small initial notch made to ensure the nucleation of a single prop- 
agating crack : no propagation, straight propagation and wavy propagation. 
Such propagating cracks were observed for velocities between 10~ 2 m.s -1 and 
1CT 5 m.s -1 . An attempt to explain this behavior was put forward [2] using 
the Cotterell and Rice criterion [6] but partially failed when the real exper- 
imental temperature profile was measured [3]. More recently the concept of 
maximum tangential stress criterion, first introduced by Erdogan et al. [7], 
was used to explain a complete set of thermal experiments [4]. Such unstable 
behavior in quasi-static crack propagation in uniaxial loading conditions was 
predicted by Adda-Bedia et al. [1] but was surprisingly never experimentally 
observed, probably because of the extreme sensitivity of this phenomenon to 
experimental conditions [1]. By working at a fixed load on DCDC specimen, 
the internal residual stresses of which were carefully relaxed by appropriate 
thermal treatment, we have shown - for the first time to our knowledge - that 
the crack front line undergoes an oscillating behavior as a result of a repro- 
ducible sequence of instabilities. The wavelength of this wavy behavior is in 
the micrometer range and its amplitude in the nanometer one. 



2 Experimental 

DCDC samples were parallelepipedic (4 x 4 x 40 mm 3 ) with a central hole of 
1mm in diameter (nominal) drilled perpendicularly to two of the large parallel 
faces. Details of experiments have been reported in a previous paper [8]. The 
sample was then placed, vertically, on its smallest face, in equilibrium on 
the horizontal basis of our experimental set-up. By means of a freely moving 
piston, the upper face of the specimen was then loaded by an increasing weight 
at a rate of about 25N. s -1 till the crack was initiated. Then the weight 
on the sample was maintained constant at a value of typically 1.5kN. These 
experiments were done at a constant temperature of 22.0 ± 1°C and at a 
relative humidity of 50 ± 5%. 

The crack front position was measured in the middle of the specimen (at equal 
distances of the lateral surfaces) by i) a standard video-recorder system for 
crack speed between 10~ 4 m.s -1 and 10~ 6 m.s -1 and ii) an optical system with 
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a magnification of x40 for lower crack speeds (10~ 6 m.s -1 - 1CT 9 m.s -1 ). In the 
latter case, the position of crack front relatively to the center of the specimen 
was measured with an accuracy of 10 fim. The stress intensity factor, Kj, as 
well as the mean crack speed, v, was deduced from the measured crack length, 
using reference [9]. The v(Ki) curve were then plotted. When the crack speed 
reached a value lower than 10~ 9 m.s -1 the samples were reloaded in order to 
gently fracture the sample all over its length. Special care was taken in order 
to preserve the two halves of the specimen from subsequent damage. 

The fracture surfaces were analyzed by an experimental system (Veeco, D3100) 
combining optical microscopy and atomic force microscopy (AFM). Both meth- 
ods could be simultaneously performed. The samples were glued on the sample 
holder. This holder was moveable under the fixed AFM-optical system using 
a stepper motor. The position of the center of the studied portion of surface 
relatively to the center of the specimen was measured with an accuracy of less 
than 10 jum. In this manner, it was possible to identify for every point of the 
AFM or optical images the value of Ki and the speed the crack had when it 
passed through this point. 

AFM experiments were done in a high amplitude resonant mode ("tapping" 
mode). Experiments were performed in ambient conditions at a relative hu- 
midity of 35 ± 5% and a temperature of 22.0 ± 0.5 °C. More details are given 
in reference [10]. 

Experiments were performed on two kinds of materials. The first one was 
soda-lime silicate glass. A thermal treatment (530°C) was done before frac- 
ture experiment in order to remove residual stresses. The second set of samples 
was made from a lithium alumino-silicate glass-ceramics. By modifying ther- 
mal treatments three types of structure were observed : the first one (sample 
A) corresponds to a pure glassy state (660° C) as identified by X-ray diffraction 
(Fig. la). The second one (sample B) is classified as a slightly unglassy state 
in which small crystals of /3-quartz phase have been nucleated (Fig. lb). This 
structure was obtained after a two-step thermal treatment at plateau tempera- 
tures of 71=750 °C and T 2 =900 °C . Samples C exhibits further crystallization 
(increased number of crystallites of larger size) due to a higher temperature 
(950 °C) for the second step of thermal treatment (Fig. lc). 



3 Results 

In the Figure 2 are plotted typical v(Kj) curves as obtained with the dif- 
ferent types of samples we studied. A noticeable feature is that for the glass 
ceramics (open symbols) the curves are shifted towards higher Kj values when 
the thermal treatment is performed. The reproducibility of these experiments 
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is revealed by the very good superimposition of the two v(Kj) curves per- 
formed on two different samples (Bl and B2: up and down triangles) which 
underwent the same thermal treatment. The fracture surfaces were then sys- 
tematically characterized by our optical-AFM system. We put a special at- 
tention to the part of the surface fracture corresponding to the low speed 
regime (10 -4 m.s -1 - 10 -9 m.s -1 ). We firstly put the stress on the lowest 
speed regime (v < 10 -8 m.s" 1 ). Scanning Electron Microsocpy (SEM) ex- 
periments revealed that corresponding fracture surfaces became increasingly 
rougher with harsher thermal treatments -from sample A (Fig. 3a) to sample 
C (Fig. 3c). These variations of RMS roughness were quantified by AFM mea- 
surements on 10 /mix 10 /im surfaces where crack speed was in the range of 
few 10 -9 m.s -1 . The results are reported on Table 1. It must be emphasized 
that RMS roughness values for pure glassy materials are very similar whatever 
the chemical composition of glasses. For the glass ceramics with an increasing 
number of larger /^-quartz nanocrystals, however, an important enhancement 
of the RMS roughness is observed. The values of RMS roughness are thus 
directly related to the length scale of structural heterogeneities. 

For higher crack speeds (except for glass-ceramics sample C; see below) an 
original phenomenon was observed (Fig. 4). Fracture surface observation re- 
vealed the presence of periodic variations in the crack path (in the direction 
normal to the mean crack plane) along the crack propagation direction. These 
corrugations in the crack path were observed : i) all along the virtual lines fol- 
lowed by both crack fronts propagating inside the specimen and symmetrically 
about the central hole, ii) on both fracture surfaces. An example (soda-lime 
glass, average crack speed v = 2.10 -5 m.s -1 and K t = 0.6MPa.m 1//2 ) is shown 
on Figure 4. The 70 /im x70 /im AFM scan of Figure 4a shows alternate black 
(deep) and white (high) bands parallel to the crack front. The sine-like shape 
of these roughness oscillations is better evidenced in Figure 4c where all the 
lines parallel to the crack propagation direction (i.e. perpendicular to the black 
and white bands) of an independent 10 //m xlO /im AFM scan (Figure 4b) 
are averaged. The occurrence (or not) of such a phenomenon was monitored 
for all the previously described samples. The following common features were 
found: 

• Samples of materials in pure glassy state (soda-lime and glass ceramics) re- 
vealed this wavy behavior for a large range of Kj values. For instance, in the 
case of soda-lime glass, the Kj domain extended between the following val- 
ues: Ki = 0.5MPa.m 1//2 and Kj = O.SMPa.m 1 / 2 respectively corresponding 
to crack front speeds of 1.5 10 -6 m.s -1 and 10 -4 m.s -1 . Between these two 
limiting values the crack ran over a distance of 2 millimeters in approxi- 
mately 5 minutes. However for samples of glass ceramics -in the pure glassy 
state- the oscillating domain is slightly narrower: AKj = O.lTMPa.m 1 / 2 . 

• One striking feature is that the A^-width of occurrence of the roughness 
waves is strongly reduced for samples in intermediate unglassy state (sam- 
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pie B): AK/ = 0.08MPa.m 1/2 . Further more no such oscillations were ob- 
served on highly unglassy samples (sample C); 

• The wavelength of these ripples, A, is in the micrometer range and its scal- 
ing with K[ values was systematically studied. Figure 5 shows the typical 
variation of A as observed in the case of soda-lime samples. It is shown that 
a large decrease of the wavelength occurs with decreasing values of Kj. Sim- 
ilar variations -in the same range of wavelength- were observed for glassy or 
slightly unglassy samples. 

• The amplitude of such corrugated oscillations is in the nanometer range. 



4 Discussion 

The experimental study presented in this paper revealed that fracture surfaces 
of various glassy materials are modulated by roughness waves the wave vector 
of which is parallel to the crack propagation direction. These roughness waves 
were observed for the various materials studied for this paper providing that 
the length scale of their heterogeneities is lower than a threshold around 10 nm. 
This length scale of heterogeneities can be very roughly evaluated by the RMS 
roughness of surface created by crack propagating at low speed (lower than few 
1CT 9 m.s -1 ), see Table 1. It means that this phenomenon is directly correlated 
to the length scale of microstructural heterogeneities of materials used for 
these experiments. 

It must be noted that a trivial AFM spurious effect can be ruled out because 
these ripples were observed too by optical microscopy for the largest wave- 
lengths (largest Kj values). It must be emphasized that, even if the aspect of 
these roughness ripples presents similarities with patterns observed in stress 
wave fractography, such an explanation is not relevant for the observations 
reported in this paper. Indeed no external transverse stress wave generator 
[11,12] was used. Furthermore the experimental set-up was working in a very 
quiet acoustical environment. Moreover the role of hypothetical residual ex- 
ternal mechanical vibrations of the DCDC apparatus may be excluded. The 
resonance frequency and the quality factor of the equivalent mechanical system 
was indeed estimated to be higher than 70Hz and 6.10 4 respectively. Conse- 
quently the hypothetical residual vibrations would have a negligible amplitude 
at the frequency range (between 0.1Hz to 1Hz) compared to the observed 
roughness waves and calculated by the ratio between the crack front speed 
and the related wavelength. 

Similar observations of wavy behavior during crack tip propagation were pre- 
viously reported [2,4]. The steady propagation of a crack localized within a 
thermal gradient was studied and an oscillatory behavior for given values 
of crack speed and glass plate width, W, was observed. The transition be- 
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tween the regimes corresponding to linear crack propagation and wavy one 
was studied as a function of specimen width for fixed temperature jump. The 
experiments by Yuse et al. [2] revealed that both the wavelength and ampli- 
tude of oscillatory behavior induced by thermal stress were in the millimeter 
range. As these length scales were of the same order of magnitude as the width 
of the region of rapid temperature change, it was assumed [2,4] that the re- 
lated instability mechanisms were only due to "external" parameters such as 
stresses induced by the thermal field and were independent of the microstruc- 
ture of broken materials. Furthermore, Yuse's experiments [2] showed that 
the related wavelength scaled linearly with W. Referring to the study Marder 
[13], it is possible to plot the variation of wavelength, A, versus Kj for the 
thermal experiments (Fig. 6). This plot reveals that log(A) slightly decreases 
with decreasing values of Ki. The data obtained in the present mechanical 
experiments can also be drawn using the same set of coordinates (see Fig. 6). 
This plot shows that a similar behavior to that revealed by thermal experi- 
ments was observed for mechanical experiments but for slightly lower values 
of Kj. This suggests that this wavy behavior could originate from the same 
phenomenon. The main difference is that, in our case, small variations of Kj 
cause larger changes in related roughness wavelength. 

A likely explanation for the occurrence of such oscillations in our DCDC ex- 
periments is that the crack front similarly undergoes a reproducible sequence 
of instabilities. Indeed such destabilizing effects in the case of uniaxial load- 
ing have already be predicted by Adda-Bedia et al. [1]. These authors argued 
that, even in quasi-static conditions, crack propagation may become unstable: 
it may deviate from straight line in a direction perpendicular to the applied 
external uniaxial loading after which a wavy crack path may appear. As al- 
ready mentioned in [1], stability thresholds drastically depend on both geo- 
metric and loading conditions. We think that our experimental set-up (DCDC 
experiments with a pure dead-weighted loading, carefully prepared samples) 
has been crucial to reveal such a phenomenon. As both wavelength and ampli- 
tude of the observed ripples are in the sub-micrometer range, the destabilizing 
factor -in the sense of Cotterell and Rice criterion [6]- from a straight prop- 
agation of the crack likely has an "internal" origin as it is connected to the 
sub-micrometric structure of materials used for these experiments. Indeed, 
when the length scale of heterogeneity is low enough such wavy propagating 
behavior is observed. On the opposite, for materials with larger heterogeneities 
the fracture surface only presents an incoherent structure (with a much higher 
roughness). The length cutoff is then in the range of lOnm. The origin of this 
instability has likely to be found in the recently predicted corrugation waves 
[14] which can be generated when fracture occurs in mode I. These authors 
indeed explain that local distortion of the crack front, due for instance to 
the presence of structural heterogeneities, may give rise to non-planar defor- 
mations of the crack front, they called "corrugation waves". As this effect 
is assumed to be due to long-ranged elastic effects causing dynamic stress 
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transfer along the crack front it likely explains why such phenomenon can be 
observed on long distances. It must be emphasized that these authors [14] 
explain that the corrugation waves may occur even in the case of a very low 
mean crack speed. Indeed, as experimentally observed [15,16] the local crack 
speed may exceed the mean crack speed by several order of magnitudes and 
lead to avalanche-like processes. Furthermore the roughness wave amplitude 
as predicted in Ref. [14] for glassy materials (8h ~ 2nm) is in the same order 
of magnitude than that observed in the current study. 

5 Conclusion 

Experiments were done on a DCDC mechanical set-up working at a fixed load. 
The internal residual stresses of specimen in different glassy materials were 
carefully relaxed by appropriate thermal treatment. In the low speed regime 
(1CT 4 — 1CT 9 m.s -1 ) we have demonstrated for the first time that the crack front 
line underwent an out-of-plane oscillating behavior along the direction of crack 
propagation. The wavelength was in the micrometer range and the peak-to- 
peak amplitude in the nanometer. This phenomenon was observed for different 
glassy materials providing that their typical length scale was lower than a 
threshold limit estimated to few nanometers. This effect was interpreted by 
the evidence of instabilities and long-range elastic effects occurring during the 
'straight' propagation of the crack front. We believed that these instabilities 
are generated by structural heterogeneities on which crack front is pinned 
and depinned. Studies are now performed at sub-micrometer scale in order to 
better understand how this reproducible sequence of instabilities is developing. 

References 

[1] Adda-Bedia, M. and Ben Amar,M., Phys. Rev. Letters (1996) 76, 1497. 
[2] Yuse, A. and Sano, M., Nature (1993) 362, 329. 

[3] Ronsin, O., Heslot, F. and Perrin, B., Phys. Rev. Letters (1995) 75, 2352. 
[4] Yang, B. and Ravi-Chandar, K., J. Median. Phys. of Solids (2001) 49, 91. 
[5] Janssen, C, In: Proc. 10th Int. Cong. On Glass (1974) pp.10.23-10.30. 
[6] Cotterell, B. and Rice, JR., Int. J. Fract. (1980) 16, 155. 

[7] Erdogan, F., Sih, G.C., Journal of Basic Engineering, ASME Transactions (1963) 
91, 764. 

[8] Despetis, F., Etienne, P. and Phalippou, J., Phys. Chem. Glasses (2000) 41, 
104. 



7 



[9] He, M.Y., Turner, M.R., Evans, A.G., Acta Met all. Mater. (1995) 43, 3453. 

[10] Marliere, C, Despetis, F., Etienne, P., Woignier, T., Dieudonne P. and 
Phalippou J., J. Non-Cryst. Solids (2001) 285, 148. 

[11] Richter, H.G. and Kerkhof, F., Fractography of Glass, Plenum Press, New York 
(1994). 

[12] Yoshida, S., Matsuoka, J. and Soga, N., J. Am. Ceram. Soc. (1999) 82, 1621. 
[13] Marder, M., Phys. Rev. (1994) E49, R51. 

[14] 1. Bouchaud, E., Bouchaud, J. P., Fisher, D. , Ramanathan, S. and Rice, J., 



J. of the Mech. and Phys. of Solids, (2002) in press - (|cond-mat/0108261 



see 



http: / /xxx.lpthe.jussieu.fr / abs / cond-mat /70108261 ) . 
[15] Delaplace, A., Schmittbuhl, J. and Maloy, K., Phys. Rev. (1999) E60, 1337. 
[16] Maloy, K. and Schmittbuhl, J., Phys. Rev. Letters (2001) 87, 105502-1. 



8 
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1 nm 


Glass ceramics : sample A 


1 nm 


Glass ceramics : sample B 


5 nm 


Glass ceramics : sample C 


10 nm 



Table 1 

RMS roughness on Wfim x lOfim fracture surfaces (mean crack speed around few 
10~ 9 m.s -1 ) for different substrates. 
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Fig. 1. X-ray diffraction patterns on glass-ceramics in a) pure glassy state, b) slightly 
unglassy state and c) more unglassy state. 



Fig. 2. Crack speed versus stress intensity factor, Ki, plots for different materials : 
Full symbols : soda-lime glass; Open symbols : Glass-ceramics (G.C.) T\ and T2 are 
the temperatures of the two-steps thermal treatment. 



Fig. 3. SEM pictures of fracture surfaces as generated by crack line propagation at 
an average mean speed in the 10~ 8 m.s -1 regime for glass-ceramics in a) pure glassy 
state, b) slightly unglassy state and c) more unglassy state. The crystalline phase 
was revealed by a mild HF etching (1%). The white marker corresponds to lOOnm. 



Fig. 4. AFM data of fracture surface of soda-lime glass; the vertical scale corresponds 
to 15 nm. a) 70/xm x 70/um scan; b) 10/um x 10/um scan; c) averaged profile of 
roughness oscillations as obtained on the Wfim x 10^m scan. 



Fig. 5. Crack speed (squares) and roughness oscillations wavelength (circles) versus 
Kj for soda-lime glass. The arrow indicates the case corresponding of AFM data of 
Figure 4. 



Fig. 6. Wavelength of corrugation oscillations versus the stress intensity factor, Kj, 
for both thermal [2] and mechanical (this paper) experiments. The dashed line is a 
guide for the eyes. 
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